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Nuclear dependence in transverse momentum distribution for Drell-Yan 
pair 

Xiaofeng Guo a * 

a Department of Physics, Columbia University; New York, NY 10027, USA 



In terms of multiple scattering picture, I compute the nuclear dependence in Drell-Yan 
transverse momentum distribution, da/dQ 2 dq T , in hadron-nucleus collisions. I present 
the results for large qx region and discuss the possible suppression in small qx region. 



1. INTRODUCTION 

Recently, it was found experimentally that the transverse momentum distribution of 
Drell-Yan pair in hadron-nucleus collisions shows nuclear dependence not linearly pro- 
portional to A, the atomic number of the target [|TJ. In addition, the average of squared 
transverse momentum, (q T ), grows approximately as A 1 / 3 [fl],0- This observed non-linear 
nuclear dependence indicates that multiple-scattering is important in hadron-nucleus col- 
lisions. In this talk, I present a calculation of double-scattering in QCD perturbation 
theory for nuclear dependence of Drell-Yan q T spectrum, da/dQ 2 dq\. 

Although double-scattering is in principle a power correction (or known as a high twist 
^ . effect) to the leading single scattering, it is not necessary small because of A 1 ! 3 enhance- 
ment from large nuclear size. In terms of QCD factorization generalized to higher twist 
0, Luo, Qiu and Sterman (LQS) developed a consistent perturbative treatment of mul- 
tiple scattering at parton level [[|. Since Drell-Yan pair does not interact strongly once 
produced, the observed non-linear nuclear dependence is a result of multiple scattering 
between the incoming beam parton and nuclear matter before the pair was produced. 
Drell-Yan da/dQ 2 dq T has two observed physical scales: Q 2 and q T . When Q 2 and q T are 
both large, da/dQ 2 dq\ is effectively an one-scale process. Method developed by LQS Q 
for double scattering can be naturally applied. However, when q T <C Q 2 , a resummation 
of large in(Q 2 /q^) is necessary for a reliable prediction of Drell-Yan q T spectrum ||. 

2. PARTON LEVEL DOUBLE SCATTERING 

Consider the Drell-Yan process in hadron-nucleus collisions, h(p') + A(j>a) — > U{q) + A, 
where the lepton pair has invariant mass (Q) and transverse momentum (qx)- As an 
example, the scattering amplitude for the lowest order double-scattering, as shown in 
Fig. |l|, has the following general form: 



M~ dx x — — F(xi,x), (1) 

J x\ — x\a + te x\ — x\b + 
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where x\ and x are parton momentum fractions, as labeled in Fig. [I], and X\ needs to 
be integrated. The sum of x\ and x is fixed by the kinematics, such as x'p f , p^ and q. 
In Eq. ([!]), the function F(xi,x) is a non-vanishing and smooth function when X\ = X\a 
and/or x\ = xib, and it is proportional to the parton fields of momentum x\p and xp 
with p = pa/ 'A. Taking the pole at x\ = x\a (or at x\ = x\b) corresponds to putting the 
propagator labeled by "A" (or "B" ) in Fig. |1] on its mass shell. In the region of interests 
(x's — 2p s ■ p > 0), both potential poles in Eq. ([I]) are in the same half of the complex 
plane. When q T ^ 0, X\a 7^ %\b- Contour integration of dx\ yields 
F{x 1A , x tot - x 1A ) F(x 1B , xtot 



M 



xib) 



M 



J- so ft— hard 



— M, 



double— hard j 



(2) 

X\A ~ XlB XlA ~ XlB 

where Xtot is the sum of the total momentum fraction from the target, and is a function 
of x', p 3 and q. Note that two amplitudes in Eq. (0) have the opposite sign. 

Amplitude M so j t _ hard corresponds to the residue at the pole "A" , at which x\ ~ and 
x ~ Xtot- We call this type of double-scattering a soft-hard scattering. In this case, the 
first scattering is effectively soft and not localized. It represents a long range correlation 
of the color field inside the nucleus. The second scattering, on the other hand, is localized 
in a distance space ~ l/x to tP ~ Because of the long range correlation of color fields, 
the soft-hard double-scattering does not have the classical double-scattering picture. Its 
contribution to the cross section can not be expressed in terms of a product of two localized 
partonic cross sections. 

Amplitude Md ou bie-hard in Eq. (g) corresponds to the residue at the pole "B". When 
qx 7^ 0, both xi = X\b and x = x to t —X\b are finite. We call this type of double-scattering 
a double-hard scattering. In this case, both partonic scatterings are hard, and they are 
localized at a distance ~ 1/xibP and ~ l/(x tot — X\b)p, respectively. Such a double-hard 
subprocess resembles the classical double-scattering picture. Its contribution to the cross 
section can be expressed in terms of a product of two localized partonic cross sections. 
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Figure 1. A sample diagram for double- Figure 2. Ratio of double-scattering to 
scattering amplitude. single-scattering. 



From Eq. (E3), we see that the double-scattering contribution to the cross section has 



three types of contributions: the soft-hard ~ M] jr t _ hard , the double-hard ~ Ml 



double— hard' 
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and the contribution from the interference terms. The interference contribution has an 
opposite sign in comparison with the other two contributions. Consequently, the inter- 
ference terms give nuclear suppression, while both the soft-hard and double-hard terms 
give the nuclear enhancement. From Eq. (0), the interference contribution is proportional 
to F*(0, x to t)F(xiB, x t ot — xib) plus its complex conjugate. When qr is large, x\b ^ 0, 
overlap in phase space for F* and F is clearly small because of the difference in parton 
momenta. However, when q? — > 0, xib ~ 0(q^/Q 2 ) — > 0. The overlap becomes large and 
the interference term becomes more important. 



3. ENHANCEMENT IN LARGE q T REGION 

In large qx region, the interference between the soft-hard and the double-hard subpro- 
cesses is not important, and consequently, the double-scattering contribution gives nuclear 
enhancement to the Drell-Yan transverse momentum distribution. The contribution from 
a double-hard subprocess can be expresses as Hj: 



da^ 
dQ 2 dq T dy 



E 

a,b,c=q,g,q 



1 



dx'f c/h (x')T» H (x a ,x b ,A) ■ — • H£ H (x b ,x') 

2x'X b S 



X 



Una 2 , 



1 



x's + u- Q 2 Q 2 



u 



)2 e « 



9Q 2 



(3) 



In Eq. (|3p, f c /h{x') is the parton distribution of parton "c" from the beam hadron h, and 
the momentum fractions of the partons "a" and "b" from the target respectively. 
The s, t, u are Mandelstam variables defined as 



(p + p'y = 2p-p'; t=(p'-q¥ 



-2p' -q + Q 2 



u 



{p-q) 2 = -2p-q + Q 2 . (4) 



In Eq. (|3|), H b D c H (x bl x') represents the calculated partonic part. Function T^ H (x a , x b , A) 
is a four parton matrix element (or correlation function) ||. For a cold nucleus, we use 
the following model : 



T° H (x a ,x b ,A) = CA^f a /N{x a )f b /N(x b ), 



(5) 



where C = 0.35/(4rQ)GeV 2 , with r = 1.1 — 1.25. At finite temperature, T^ H depends 
on the density matrix at the given temperature. 

For a soft-hard subprocess, the contribution has the following general form||: 



da^ 
dQ 2 dq T dy 

with 



E Jdx'f c/h (x')^ s ac H (x a ,x',A).^ 



12na 2 



s x's + u — Q 2 



12 e a » 



9Q 2 



(6) 



d 2 ( 1 SH 



dx\ \x a 



2q 2 



{x's + u-Q 2 ) 2 



+ 







dx a \ x a 



2(Q 2 



u 



x's(x's + u — Q 2 ) 



(7) 
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In Eq. (|7|), is the calculated partonic part, and T a (x a ,A) is the soft-hard matrix 

element. For a cold nucleus, LQS proposed the following model |4]]: 

T a (x a ,A) = X 2 A^f a/N (x a ), (8) 

where A 2 ~ 0.05 — 0.1 GeV 2 , and A is estimated from the data on di-jet momentum 
imbalance |J. 

Combining contributions from double-hard and soft-hard subprocesses, and using the 
models for T^ H and T a , large nuclear enhancement for Drell-Yan production in hadron- 
nucleus collisions was obtained ||. A typical enhancement in large region is shown in 
Fig- 1- 

4. SUPPRESSION IN SMALL q T REGION 

When q T — > 0, x\a ~ X\b ~ 0, both soft-hard and double-hard contributions become 
divergent. But these divergences are canceled by the interference terms from Eq. (fj), 
because 

M (F{x XA ^-x XA )\ _{F(x 1B ,x t ^ t -x 1B )\ ^ fin . te; (g) 

\ %IA %1B J soft-hard \ XlA %lB /double-hard 

as qT —> 0. In addition, similar to the single-scattering, double-scattering subprocesses 
develop the collinear and soft divergences when qT — > 0. Although all divergences are 
canceled after including necessary virtual diagrams and proper collinear subtractions, 
there are large £n(Q 2 / q T ) for every power of a s . A systematic resummation at the presence 
of double-scattering is necessary for a reliable prediction in small g-r region ||. 

Because of very small A-dependence in Drell-Yan da/dQ 2 , and large nuclear enhance- 
ment for large qT region, one expects nuclear suppression in small qT region. This is consis- 
tent with above discussion on the role of quantum interference between the soft-hard and 
double-hard subprocesses. It is known that nuclear shadowing in the parton distribution 
of momentum fraction x is due to quantum interference of multi-parton recombination 
in the longitudinal direction 0. If the same thought is applied to the transverse di- 
rection, a nuclear suppression in small qT region is expected from quantum interference 
between soft-hard and double-hard subprocesses. The expected ratio of double-scattering 
to single-scattering contribution in small qT region is sketched in Fig. | [[|. 
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